Association between homocysteine and coronary artery disease—trend over time and across the regions: a systematic review and meta-analysis

Background The association of homocysteine with coronary artery disease (CAD) has been explored previously with mixed findings. The present Systematic Review and Meta-Analysis (SRMA) has assessed the pooled estimate of association between homocysteine (Hcy) and CAD, and its variation over the period and geography. Methods Systematic literature search was done in PubMed, Scopus and Cochrane to identify the observational studies that have reported mean Hcy among cases (CAD) and control. The SRMA was registered in PROSPERO (ID-CRD42023387675). Results Pooled standardized mean difference (SMD) of Hcy levels between the cases and controls was 0.73 (95% CI 0.55–0.91) from 59 studies. Heterogeneity was high (I2 94%). The highest SMD was found among the Asian studies (0.85 [95% CI 0.60–1.10]), while the European studies reported the lowest SMD between the cases and controls (0.32 [95% CI 0.18–0.46]). Meta-regression revealed that the strength of association was increasing over the years (Beta = 0.0227, p = 0.048). Conclusions Higher homocysteine levels might have a significant association with coronary artery diseases, but the certainty of evidence was rated low, owing to the observational nature of the studies, high heterogeneity, and publication bias. Within the population groups, Asian and African populations showed a greater strength of association than their European and American counterparts, and it also increased over the years. Supplementary Information The online version contains supplementary material available at 10.1186/s43044-024-00460-y.


Background
Cardiovascular disease (CVD) is a major public health and clinical problem across the world.Though mortality due to CVD is decreasing in developed countries, the proportion of deaths due to coronary heart disease (CHD) is increasing (of all the deaths among people over 35 years of age, around one-third are due to CHD) [1].Over ninety percent of CAD events occur in individuals with at least one risk factor.Among the risk factors, few are non-modifiable, like increasing age, male gender, history of premature CAD among first-degree family members, and some genetic factors, while many are modifiable or preventable risk factors like hypertension, high fasting plasma blood glucose, obesity or overweight, physical inactivity, dyslipidemia (high low-density lipoprotein and low high-density lipoprotein) cigarette smoking, stress, suboptimal diet, etc. [2] But in recent years, studies have reported that known classical risk factors can explain only half to one-third of atherosclerotic vascular events [3].Certain novel potential risk factors that are being studied and hypothesized are high-sensitivity C-reactive protein, lipoprotein (a), plasma fibrinogen, plasma homocysteine, plasminogen inhibitor type I, endogenous tissue plasminogen activator (tPA), estrogen deficiency, microalbuminuria, increased levels of the leukocyte enzyme myeloperoxidase, etc. [2,3] Homocysteine (Hcy) has been studied as an independent risk factor for vascular disease.Researchers across the world are interested in finding out its association with CAD, cerebrovascular disease, and peripheral vascular disease.Fasting plasma Hcy concentrations between 5 and 15 µmol/L are considered normal [3].The correlation between hyperhomocysteinemia and atherosclerotic disease was first detected more than 50 years ago by McCully in 1969 [4].One of the earlier meta-analysis done by Boushey et al. concluded that a 5 µmol/L increment in total homocysteine level and 0.5 µmol/L (20 mg/ dL) increase in serum cholesterol elevate the similar risk of CAD [5].Arnesen et al. observed a positive relationship between homocysteine and the risk of myocardial infarction in their large prospective study [6].Schnyder G. et al. suggested that total homocysteine, along with age and gender, strongly predicts the severity of coronary artery disease and should be assessed for the CVD risk profile of patients as an independent cardiovascular marker [7].On the contrary, a few studies did not find a significant association between homocysteine and coronary artery disease [8][9][10].Thus, the studies have reported mixed findings in terms of the relationship between plasma homocysteine and cardiovascular disease.Preliminary systematic search showed that previously published SRMAs in 2008 [11] and 2022 [12] reported a significant association between Hcy and CHD.
However, these analyses had the following limitations: The SRMA published in 2008 included studies representing population from North America and European regions only, undertook Medline and Cochrane database search only [11], while three primary database search is recommended for SRMAs [13].The outcomes of the included studies were also diverse like any CHD event including CHD death, MI, revascularization procedures, CVD death or stroke [11].In the 2022 SRMA, the search strategies (database specific), inclusion-exclusion criteria and the outcomes were not comprehensively defined and/ or published.Publication bias was not assessed quantitatively, and outlier determination was not done [12].Certainty of the evidence was not determined [11,12].Since studies considered different cut-off values of Hcy based on median, tertiles, quartiles or quintiles, risk ratio (RR) was estimated for each study based on assumed log-linear association between CHD risk and Hcy.So, actual RR of the studies was not used in meta-analysis.Based on the above critical gaps identified, we undertook the index meta-analysis to determine the association between Hcy and the CAD by adopting the following comprehensive and transparent methodology.

Systematic search
The index SRMA adhered to the PRISMA guidelines (Additional file 1: Table S1).The research question for the index SRMA was: What is the association between the plasma homocysteine level and the occurrence of coronary artery disease?which is elaborated in Additional file 1: Table S2.An extensive literature search was done in various databases: PubMed, Scopus and Cochrane.The authors also searched Google Scholar to find other related articles.The base search strategy was formulated for PubMed.The same search strategy was used for the rest of the database as per the required format (Additional file 1: Table S3).We also reviewed the references of the eligible articles and found more studies that could be included in this systematic review and meta-analysis.

Selection of study
After the removal of duplicates, screening of the studies was done based on the criteria elaborated in Additional file 1: Table S2.Titles and abstracts were screened for the outcome of interest, i.e., Acute Coronary Syndrome or Coronary Artery Disease, and the exposure of interest, i.e., total plasma homocysteine concentration.Through title and abstract screening, observational studies like cross-sectional, case-control and nested case-control studies were identified where the case group represents acute coronary syndrome (ACS) or CAD patients and the control group comprises healthy individuals or patients free from coronary artery disease and exposure, i.e., plasma homocysteine level, has been measured.This initial title-abstract screening of the studies was done independently by two reviewers (SVU and AG).Full-text articles of these short-listed studies were retrieved and again screened by the same two authors (SVU and AG) against predefined study eligibility criteria.At any level, discrepancies about the eligibility of the study were discussed and resolved by the reviewers, and adjudication was sought from the third reviewer (BKP), if consensus was not achieved between the two reviewers.

Data extraction
All the relevant data were extracted independently by the two authors (SVU and AG).For that, a common data extraction format was prepared using Microsoft Excel.Information like the name of the author, publication year, study type, number of cases and controls, and characteristics of cases and controls were retrieved.Information about exposure, i.e., plasma homocysteine, was collected, like mean and SD among case and control groups, fasting or non-fasting blood samples and methods of homocysteine estimation.If the data were reported as a mean and standard error in the study, then SD was calculated.Other relevant information was also recorded, like the mean age, age group and gender of the study participants.

Assessment of the quality of the studies
To assess the study quality, the National Heart, Lung, and Blood Institute (NHLBI) tool was used, which was developed for observational studies [14].Two reviewers assessed the quality of the studies and third reviewer resolved if any discrepancy or disagreement raised.Studies fulfilling at least 75% criteria were labeled good-quality studies, while those with 50% to 75% were fair, and less than 50% were considered poor-quality studies.

Data analysis
The pooled estimates of the outcomes, along with a 95% confidence interval (95% CI), were measured in terms of standardized mean difference (SMD).I 2 statistics were applied to assess the heterogeneity of studies, and I 2 > 50% was considered substantial to high heterogeneity [15].If I 2 > 50% was found, a random-effects model (the Der Simonian and Laird method) [16] was applied.The prediction interval of the pooled estimate was determined based on the Tau 2 statistics [17].Heterogeneity was explored by undertaking subgroup analyses based on the following variables: case group (ACS/CAD), geography (continent of study), gender (males and females/ males only), type of blood sample (fasting or non-fasting), and study period.Assessment of publication bias was planned by means of the Funnel plot, Doi plot, and LFK index if more than ten studies were found eligible for meta-analysis.Baujat plot and influential plots were done to identify the influential studies.A leave-one-out analysis was planned to estimate the impact of each study on the pooled outcome estimate and the heterogeneity.A sensitivity analysis was planned after eliminating the poor-quality studies and influential studies.A p-value of < 0.05 was deemed significant.

Certainty in the evidence
The study evaluated and summarized the pooled estimate's certainty for each outcome using the Grading of Recommendations, Assessment, Development, and Evaluations (GRADE) methodology [19].
The systematic review and meta-analysis was registered with PROSPERO (ID-CRD42023387675).

Eligible studies
The databases were thoroughly searched, resulting in a collection of 2504 distinct records.These records were then subjected to title-abstract screening, leading to the identification of 379 studies that were deemed suitable for a comprehensive review.After conducting a thorough examination of the full texts of these studies, 59 of them were determined to be eligible for data extraction and were thus included in the systematic review and metaanalysis (Fig. 1).

Risk of bias assessment
Upon assessing the risk of bias, it was determined that out of the 59 studies included in the meta-analysis, 55 were classified as being of good or fair quality, while four studies were of poor quality (Additional file 1: Table S4).

Outcomes
The pooled SMD of the homocysteine levels between the cases and controls was 0.73 (95% CI 0.55-0.91),which implies a significantly higher homocysteine levels among the cases than the controls.There was high heterogeneity between the studies (I 2 = 94%, p < 0.01).The prediction interval was between − 0.60 and 2.06 (Fig. 2).Sensitivity analysis following the removal of poor-quality studies resulted in a similar pooled SMD of 0.73 (95% CI 0.54-0.91)(Additional file 2: Figure S1).
Two of the 59 studies were identified as potential outliers (Additional file 2: Figure S2), and the sensitivity analysis revealed a slightly lower SMD of 0.64 (95% CI 0.50-0.78).Although there was a minor decrease in heterogeneity, it still remained at a high level with an I 2 value of 92% (p < 0.01).Additionally, the prediction interval became narrower, ranging from − 0.35 to 1.63 (Additional file 2: Figure S3).
The leave-one-out analysis revealed that there was not much variation observed in the heterogeneity, with I 2 values ranging from 93 to 94%.Similarly, the pooled standardized mean difference (SMD) estimate remained relatively stable, varying from 0.68 to 0.75 throughout the analysis (Additional file 2: Figure S4).
Table 2 presents the pooled estimate resulting from the subgroup analysis.Except for the Europe subgroup (I 2 51%), heterogeneity remained substantial and even among all the subgroup analysis.Diagnosis (ACS/CAD) and the fasting status of the blood sample did not significantly affect the SMD between the cases and controls.There was a significant difference in the homocysteine levels association with the CAD/ACS between different regions of the world, with the highest SMD among the Asian studies (0.85 [95% CI 0.60-1.10]),while the European studies reported the lowest SMD between the cases and controls (0.32 [95% CI 0.18-0.46]).Similarly, the  time period of the study had a significant impact on the pooled SMD of the Hcy levels, with the post-2000 SMDs being significantly higher than the pre-2001 levels.Metaregression also revealed that the strength of association was increasing over the period in the index meta-analysis (Beta = 0.0227, p = 0.048).(Fig.3).Relatively, studies that included only males had a significantly lower SMD in the homocysteine levels (0.40 [95% CI 0.28;0.53])than the    2).Funnel plot and Doi plot asymmetry analysis indicated the presence of potential publication bias.This bias was statistically represented by the LFK index, which yielded a value of 2.69 (Fig. 4a, b).
The certainty of the evidence was rated low, owing to the observational nature of the studies, high heterogeneity, and publication bias (Table 3).

Discussion
Of the novel biomarkers of CAD, plasma homocysteine is being studied to determine whether it is an independent modifiable risk factor or not.The mechanism by which high concentrations of plasma homocysteine lead to atherothrombosis is endothelial cell damage impairing its function, which has been shown in many in vitro studies as well.High Hcy levels increase oxidative stress, inducing inflammatory processes, increasing endoplasmic reticulum stress and apoptosis, increasing autoimmune reactions, and ultimately increasing the coagulation cascade [68,78].The index meta-analysis, which is the first to present the pooled estimate of this association, included a total of 9381 cases of CAD and 12,188 controls from the 59 studies and revealed a significantly higher plasma homocysteine concentration among the cases (SMD of 0.73, 95% CI: 0.55-0.91).This suggests a significant association of homocysteine with the occurrence of CAD, although the overall certainty of the findings is low, mostly due to high heterogeneity.The high heterogeneity in the association reported in the index meta-analysis could be due to the difference in inclusion and exclusion criteria of cases and controls, inadequate sample size, different methods of homocysteine measurement and fasting status of blood sample collection, the presence of other traditional risk factors, and the wide range of time periods of the included studies.Plasma homocysteine concentration is affected by many factors like age, gender, ethnicity, nutritional deficiencies like folate and Vitamin B12, and renal and liver function [79].In addition to these factors, the enzyme methylene tetrahydrofolate reductase, encoded by the MTHFR gene, also regulates the homocysteine level.This gene metabolizes and removes homocysteine by using folate.A polymorphism of the gene MTHFR C677T reduces the efficiency of the enzyme and leads to an increase in plasma Hcy concentration [80], which might have also contributed to the heterogeneity.Twenty-one of the 59 studies have investigated the different genetic and environmental interactions involved in the occurrence of   [69].Studies adopted different case groups, like patients with acute conditions or those who had coronary events in the past, where the diagnosis was based on the extent of stenosis on angiography.The pooled estimate of SMD was high among the CAD group, where patients with any degree of stenosis were included.Heterogeneity was little reduced (from 95% in the CAD group with ≥50% stenosis, 94% in the ACS group and 92% in the CADcriteria unspecified group).Various studies reported that patients with severe CAD in terms of the number of vessels affected or higher Gensini scores had significantly elevated mean total homocysteine levels [23,82].But Bozkurt A et al. reported that homocysteine concentration was unrelated to the extent (in terms of the number of vessels affected) and severity of the disease [34].Pooled SMD is high for the studies done in Asian countries (pooled SMD 0.85, 95% CI: 0.60-1.10),followed by African countries (pooled SMD 0.75, 95% CI: 0.18-1.32).Heterogeneity is also high among studies done in countries on these two continents, while studies in European or American regions showed less heterogeneity as well as a small pooled SMD.Since homocysteine is believed to be determined by nutritional deficiencies like folate and Vitamin B12 as well as ethnicity, the finding of the index subgroup analysis can be explained by studies exploring the folate levels in various geographical regions.A meta-analysis of MTHFR polymorphism with CHD risk done by Clarke R et al. reported folate levels as low in Asian and European un-supplemented populations, intermediate folate levels in the supplemented European population and un-supplemented US and Australian populations, while high folate levels in the supplemented US and Australian populations [80].Chambers JC et al., in their two parallel case-control studies done among two ethnic groups, European and Indian Asian men residing in Europe for an average of 27 years, reported that fasting Hcy concentrations were high in Indian Asian men compared to European men, and the age-adjusted difference was 6% [83].Another study conducted to understand the impact of migration on the risk of coronary heart disease revealed that serum homocysteine levels were significantly higher among the Indian participants residing in India compared to Indian-origin participants who migrated to Sandwell, UK.This was consistent with low serum folate levels among Indian participants compared to migrated participants [84].CAD is one of the major public health concerns across the world and is known for its multifactorial causation.Thus, genetic as well as environmental factors (folate levels) have an interactive effect on the Hcy levels.
Subgroup analysis based on gender showed that pooled SMD was lower in studies that included only male participants compared to studies with both genders.Studies have shown that males had significantly higher Hcy values than females at each age range, which could be a contributing factor to gender differences in developing CAD [34,85,86].The majority of the studies have documented fasting homocysteine estimation through high-performance liquid chromatography, fluorescence polarization immunoassay, or ELISA methods.Subgroup analysis shows a lower but insignificant pooled SMD when homocysteine is estimated from a non-fasting blood sample.The estimated pooled SMD was lower in the studies published before 2001 (0.37, 95% CI: 0.23-0.52),with substantial heterogeneity.Almost all the studies were done in Europe or the American region except one, which was from Japan.After 2010, all the studies were published in Asian or African countries.So, this difference in pooled estimates over the period of time might be due to the region of study.It is also possible that over the span of three decades, lifestyle has changed considerably and all the known risk factors have become more prevalent, thus increasing the vulnerability for coronary artery disease.Food habits and probably the quality of food have also changed over time, so nutritional deficiencies particularly folate and Vitamin B12 deficiencies could have increased, which in turn is associated with a rise in homocysteine levels, particularly in low-and middle-income countries in Asia and Africa.

Strengths and limitations
The index meta-analysis is the first study to present the global, regional, and temporal pooled association estimates of the Hcy with CAD.The quality of the studies was evaluated by standard tools, and sensitivity analysis was undertaken to improve the robustness of the findings.Heterogeneity was explored by means of appropriate and feasible subgroup analyses.The GRADE profile was applied to present the certainty of the evidence from the meta-analysis.However, the index meta-analysis had the following limitations: High heterogeneity between the studies persisted even after subgroup analysis might be due to the presence of publication bias, language bias or genetic reasons.Although we searched three of the major databases (PubMed, Scopus and Cochrane) along with Google Scholar, more databases could not be included, because of limited resource availability.

Conclusions
Overall, even though higher homocysteine levels might have a significant association with coronary artery diseases, the certainty of evidence is low.Within the population groups, Asian and African populations showed a greater strength of association than their European and American counterparts.High heterogeneity, which is a major factor impacting the certainty of the evidence, needs to be explored in future studies by reporting and conducting subgroup analysis based on determinants such as genetics, sex and folate levels.Primary studies can be designed to alleviate all probable confounding factors to assess the predictive role of homocysteine in CAD and whether it is a modifiable risk factor.

Fig. 1
Fig. 1 PRISMA flowchart of the study

Fig. 2
Fig. 2 Forest plot showing the pooled estimate of the standard mean difference in the homocysteine levels between CAD and controls

Fig. 3
Fig. 3 Meta-regression of the outcome measure over the time

Fig. 4 a
Fig. 4 a Publication bias assessed by funnel plot and b publication bias assessed by Doi Plot

Table 1
Characteristics of the Individual studies included in the meta-analysis (n Author (

Table 2
Pooled estimates of Hcy SMD according to the subgroup analysis

Table 3
Certainty of evidence on the association between the homocysteine and CAD by GRADE profile